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REMARKS  ON  THE  REFLECTION  OF  CONICAL  SHOCKS 
By  K.  0.  Friedrichs 

The  theory  of  reflection  of  e  plane  shock  at  a  plane  wall  has  been 
thoroughly  treated  by  v.  Neumann.*  Only  for  sufficiently  small  angles  of 
incidence  CK   (see  Fig.  1)  is  s  "regular"  reflection  possible  while  for 
larger  angles  CX.  up  to  90°  a  Mach  reflection  will  occur  (see  Fig.  2). 
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Fig.  1  Fig.  2 

It  is  naturally  of  interest  to  investigate  whether  or  not  the 
situetion  is  similar  for  c.xiclly  symmetric  shocks.  Regular  reflection 
at  an  axis  would  involve  two  conical  shocks,  the  incident  and  the  reflected 
one,  both  having  the  tip  in  common  (see  Fig.  3).   The  flow  is  assumed  to 
be  exially  symmetric  and  steady  when  observed  from  the  tip.   It  is  natural 
to  make  in  addition  the  following  basic  assumption;   along  each  straight 
ray  through  the  tip,  pressure,  density,  and  even  velocity  are  constant. 
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Fig.  3 
It  is  the  aim  of  the  present  remarks  to  show  that  under  this  assumption 
no  regular  reflection  of  conical  shocks  is  possible.   It  is  an  open  question 
whether  Mach  reflection  takes  place  instead  of  regular  reflection  even  for 
small  angles  of  incidence,  or  whether  the  basic  assumption  should  be 
abandoned.   Experimental  evidence  of  regular  conical  reflection  seems  to 


*  A  preliminary  draft  of  v,  Neumann's  report  was  available. 
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point  to  the  second  alternative, 

Axially  symmetric  flow  md  shocks  satisfying  the  basic  assumption  have 
been  treated  by  Buse-n^nn,  by  Taylor  end  Maoooll  and  by  others.*  The  main 
objective  of  those  authors  was  to  investigate  the  flow  which  results  when 
a  conical  projectile  moves  with  super-sonic  speed.   The  nroblem  as  a  whole 
is  explained  in  a  most  lucid  wcy  by  Taylor  and  Maccoll.  As  to  the  analytical 
formulation  of  the  problem  Busemann's  procedure  is  more  elegant  and  flexible. 
We  shall  follow  Busemr.nn's  procedure  in  the  present  memorandum. 

We  first  derive  shortly  tl^e  differential  equation  for  the  axial ly 
symmetric  flow  assuming  that  pressure  p  and  density  p  are  connected 
through  the  adiabatio  relation 

(i)  p  -  k/>r,  y  - 1.4. 

Let  x  and  r  be  cylindrical  coordinates  and  let  u  and  v  be  axial 
and  radial  components  of  the  velocity.   The  irrot  tional  character  of  the 
flow  is  expressed  through 

(2)  vx  =  ur  , 

while  the  continuity  equation  is 

(3)  r(/Ju)x    (rPv)r   *   0. 
By  virtue  of  Bernoulli's  equation 

(4)  Hu2  +•  v2)   -t-   _>T„.  p//)    i?  °  2   "   const. 

and  relation  (l),  one  may  eliminate  the  density  jO    from  equation  (3)  and 


*  G.  J.  Taylor  and  J.  W,  Maccoll.   The  air  pressure  on  a  cone  moving  at 
high  speeds,  proc.  R.  Sec.  A.  Vol.  139  (i933)  t>.  278.     (| 
A.  Busemann,   Die  Achsensymmetrische  Kegelige  Uborschallstromung. 
Luftfahrtforschung,  Vol.  19  (1942)  p.  137,  and  literature  given  there, 
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obt?..ins 


(5)  (1   -  «   )   ux  +  (1   -  r.)   vr  +    v     .  2  ^  v     ~-    O       . 

c  c6 

where 


(6)  c'     a    yyv>//> 

is  the  velocity  of  sound.     By   (4)   we   have 

(7)  C/2=..21.-.i(^.u2-v2). 

The  basic  assumption  implies  that  u  end  v  depend  only  on 

(8)  t  =  x/r  , 
Equation  (2)  becomes 

(9)  vt  -r  tut   =  0 

while  (5)  reduces  to 

2  2 

(10)  (1  -  "   )  ut  -  (1  -  I_  )  t  vt  +•  v  -  2  ™  vt  =0. 

c4  c^  c- 

The  variable  t  can  be  eliminated  by  (9)-  it  is  further  possible  to  elim- 
inate differentiation  with  respect  to  t   and  to  consider  v  as  a  function 
of  u  .   To  this  end  we  write  equation  (9)  in  the  form 

(11)  t   =  -  vt/ut  =  -  vu  j 
Differentiation  of  (11)  with  respect  to  t  ^ives 

(12)  1       »  "vuuut   ' 

Expressing  u^  ,   v,  ,  t  through  v   and  vyu  ,  we  obtain  from  (10)  ' 

(13)  (1  -  _u5  +  (1  -  4^  vu  -  v  vuu  -  2  ^  vu  =  °. 
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which  can  also  be  written  in  the    r'orm 


(14) 


w        =•    1 
uu 


(u    +•  vru)s 
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This  is  Busemann's  equation*  (derived  in  a  different  manner). 

Jt  is  dear  from  (9)  and  (10)  that  u^  ^   0,  for  every  solution  of 
(8)  end  (10)  as  long  as  u^O  ,  (what  happens  at  u  =0  requires  special 
investigation),  Every  such  solution  thus  leads  to  a  solution  of  (14)  with 


(15) 


uu 


?t  0  . 


On  the  other  hend,  every  solution  of  (14)  which  satisfies  (15)  leads 
to  :..  solution  of  (9)  ,nd  (10)  when  t  is  defined  by  (11).   In  other  words, 
every  section  of  a  solution  of  (14)  for  which  vu  is  monotone  leads  to  a 
possible  flow  when  the  ray  to  which  u  and  v  are  to  be  attached  is  given 
by  (11)  and  (8). 


*  It  may  be  mentioned  that  Busemann  gives  an  elegant  geometric  interpretation 
of  equation  (14) : 

(16)  R  -  N/(l  -  .Hi)  . 

c 
R  being  the  radius  of  curvature  end  the  meaning  of  N  .and  U  being  obvious 
from  the  figure . 
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When  v  as  a  function  of  u  is  represented  by  :  curve  in  the 
(u,v)-plane,  the  direction  of  the  ray  assigned  to  a  point  (u,v)  is  that 
of  the  normal  to  the  curve  at  this  point;  as  seen  from  (11). 

The  relations  governing  the  transition  through  a  conical  shock  arc 
the  some  as  for  the  plane  obliaue  shock;  the  curvature  of  the  shock  cone 
does  not  enter.  When  the  shock  cone  is  a  strr ight  cone,  as  is  assumed, 
the  jumps  of  u,v  of  p,  ,   and  of  the  entropy  are  constant  along  each 
ray  when  the  basic  assumption  is  setisfied  on  one  side;  consequently  this 
assumption  remc ins  satisfied  on  the  other  side.   The  flow  may  continue 
as  a  pctentiel  flow  with  constant  entropy  after  crossing  the  shock.   In 
other  words  the  assumption  of  proper  conical  shocks  is  compatible  with 
the  basic  assumption. 

Suppose  a  flow,  characterized  by  p  ,  /)  ,  u  ,  v   crosses  such  a 

conical  shock,   (it  is  to  be  noted  that  this  can  occur  only  if  the 

velocity  w   =  xfi    +•  v?-  is  supersonic,  i.e.  if  wQ  >  eQ) .   The  velocity 

(u,  ,v,)  immediately  past  the  siiock  is  located  on  the  loop  of  a  strophoid 

in  the  (u,v)-plane.  The  inclination  of  the  ray  which  generates  the  shock 

cone  is  perpendicular  to  the  straight  connection  between  (u0  ,  v0)  and 

(u-,v, ).   T^e  positions  of  the  cones  corresponding  to  the  cases  (i); 

v.,  >  v   and  (il):  v  <  v   are  indicated  in  Figure  6. 
1   o  1   o 
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Case  (I)  Case  (II; 

Fig.  6 
If  the  state  (l)  past  the  flor  is  prescribed,  the  points  (uo,vQ) 
corresponding-  to  possible  states  In  front  of  the  shock  *re  situated  on 
the  "tail"  of  the  strophoid  through  f^.Vj)  as  indicated  in  Figure  7  for 


case  (I)  . 
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PigV  V 
When  the  flow  on  either  the  front  or  the  back  side  of  the  shock  is 
to  be  continued  according  to  differential  equation  (14),  the  slope  of  the 
(u,v;-curve  is  to  be  so  determined  that  the  ray  given  by  (11)  coincides 
with  the  shock.   Since  this  rsy  is  bo  be  normal  to  the  (u,v) -curve  on  the 
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one  hand  and  perpendicular  to  the  straight  segment  connecting  (u  ,v  )  with 
(u  ,v  )  on  the  other  hand,  the  (u,v)-curvo  should  begin  or  enter  in  the 
direction  of  this  segment.  The  slope  of  the  (u,v) -curve  is  thus  given  by 

(17)      vu     =     ;v:   -  v0)/(Ul  -  uj  . 

The  discussion  of  conical  shocks  by  Busemenn  end  by  Taylor  and  Macooll 

was  restricted  to  case  (I)  with  u  =  w  >  0  and  v   «  0  .   This  case 

0    0  o 

occurs  when  a  constant  axial  flow  is  deflected  by  a  conical  projectile. 
It  may  shortly  be  indicated  how  Buseirurn  treats  this  problem.  Through 
tve  shack  transition  relations  the  flow  velocity  (u  ,v  )  past  the  shock 
is  givenj  (observe  that  the  third  transition  relation  guarantees  that  the 
Bernoulli  constant  §•  ac     is  the  same  before  and  after  the  shock).  A 
solution  of  equation  (14)  is  to  be  found  whose  c-raph  passes  through  the 
point  (u,  ,v.).  The  slope  vr  of  this  curve  is  given  by  (17).  The  solution 
is  no":  to  be  so  continued  that  t  =  x/r  increases,,  i.e.,  in  view  of  (11), 
v   decreases  up  to  a  point  at  which  the  flo,r  and  the  ray  have  the  same 
direction,  i.e.,  where  v/u  —  x/r,  or  where  the  normal  passes  through  the 
origin;  such  a  point  may  be  called  an  end  point.  This  end  point  depends 
on  the  choice  of  the  point  (u. ,v. )  on  the  strophoid.  The  manifold  of 
endpoints  that  can  be  reached  from  (wo,0)  forms  a  curve  which  Busemann 
calls  "ao'-'le  curve"  in  view  of  its  peculiar  shape,  see  Fig.  8. 
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Ir  this  procedure  the  shock  is  nrescrifce-'  and  the  end  direction  is  found. 
If  the  end  direction  is  prescribed  one  may  find  the  corresponding  point  on 

apole  curve  by  intersecting  it  with  the  appropriate  ray  through  the 
origin.  In  general  there  will  be  two  intersections  of  which  the  one  with 
the  weaker  shock  is  likely  to  occur  in  reelity  parallel  to  projectile.  * 
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We  nov.  return  tc  our  problem  of  regular  conical  reflection.   The 
incident  shock  is  of  the  first  type,  the  state  (0)  in  front  of  it  being 
constant  with  axial  flow  direction,  u  =  v    >  0  ,  vQ  -  0  .  The  flow 
past  the  incident  shock  (l)  will  point  toward  the  axis,  v,  <  0  ;  it  will 
not  remain  const- nt  and  hence  the  state  (?)  in  front  of  the  reflected 
shock  .vil1  differ  from  state  (l).   The  reflected  shock  is  of  type  (il); 
the  flow  behind  it  (3)  may  cither  be  immediately  axial  and  constant  or** 
it  will  vary  until  its  direction  has  become  axial  and  be  constant  from 
there  on  (4),  of.  Figure  10. 


It  may  bo  mentioned  that  in  the  procedure  of  Taylor  and  Maccoll  one 
is  to  begin  with  the  end  direction:  the  shoe1  is  then  found  by  foil 
ing  the  solution  of  (10)  backwards.  This  procedure  has  advantages 
-  icn  single  cases  are  to  bo  investigated. 


OY 


*•'••   This   possibility  was   pointed   out  to  me  by  v.   Neumann. 
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Fig.  10 

The  flow  between  the  two  shoots,  i.e.  between  (l)  and  (2),  (>.nd 

also  that  between  (3)  ;n'i  (4)  )  is  ."warned  by  the  difforentirl  equatii 

(14)  for  v  as  c.  function  of  u  .  The  initial  direction  at  (l)   (and 

also  at  (3)  )  is  determined  from  relation  (17).   In  order  to  establish 

the  impossibility  of  suoh  ?.  co     r  tion,  we  shell  show  that  the  expression 

Tr  v    is  positive  pent  any  shock,  hence  at  (l),  while  it  is  negative  in 

front  of  any  shock,  hence  rt   (2).      The  sirn  of  v  vuu  must,  therefore, 

change  in  the  transition  from  (1)  to  (9).  A  change  of  the  sign  of  vuu 

is  forbidden  by  (15).  That  v  also  cannot  change  sign  can  be  seen  as 

follows;  In  state  (l)  we  have  v  <  0  ,  hence  vuu  <  0  .  If  vuu  does 

not  change  sign,  v  u  <  0  when  v  "0  is  reached.  Clearly,  vy  ^  0 

when  v  =  0  is  reached  for  the  first  time.  On  the  other  h_.nd,  v  =■- 0, 

v  i-   0  implies  v   ?-  0  ;  this  is  seen  by  differentiating  (14)  and  then 
u  uu  ' 

setting  v  ■•"  0  *,  Thus  a  contradiction  is  established. 

We  prove  our  statement  concerning  v  vuu  by  direct  calculation. 
Let  (u  ,v  )  correspond  to  any  state  which  may  b'  in  front  of  or  behind 
a  shock;  see  Figure  11.  Let  (u,v)   correspond  to  the  state  on  the  other 


*  The  justification  of  this  argument  follows  from  the  expansion  of  v(u) 
in  the  neighborhood  of  v  =  0  given  by  Busemann,  cf.  his  (14). 
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side  of  the  shock.  Let  (u  ,v  )  be  the  coordinates  of  the  point  (u,v)  with 
respect  to  a  system  which  is  obtained  by  rotating  the  u-axis  such  that  it 
passes  through  the  point  (uo»v0) •   Then  we  have  (cf .  Buscmann) 


(13)    vl      =   (w0   -  u/ 


w. 


W0 


with  the  notations 

(19)    wo  =  yu*  ^v2o , 

^2  / 
Wj   =  1/  c  /w0  , 

w2   x  wx  + 

y  +  i  ' 

/  t  1 

X 


w 


The  slope  v   of  the  solution  of  (14)  which  begins  or  ends  with  the 
state  (0)  is  given  by  (17).  Hence  ve  hr vc  from  (18) 


(20) 


w2 


vr, 
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is  eesilv  determined  to  be 


U2  =      K  -  «*)2    w2 

or  from  (18) 

(21)  U2  -   ^2  -  u»)  wo 

X 

In  view  of  (7)  and  (19)  we  have 

(22)  Ao   =  wQ  (wj  -  x/  wQ  )  , 

Hence,    from  (21)    and   (22), 

II2  2 

1  -     -?      -      -  K  - «*)    I°_    > 

c  \£ 

or  by  virtue   of  (14) 

(23)  v  vm,      =     -  (wo  -  u$)  *§  (H-^  )      ' 

"X  c^ 

Now,    if  the   state   (0)    is   in  front  of  the   shock  we  have     u„.  <.  w     ,   hence 

'  *    o 

v  v      <  0;  if  (0)  is  past  the  s"-ock  we  have  u   >  w   .  hence  v  v   >  0, 
uu  *     0  uu 

Thus  our  statement  is  proved. 
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